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ABSTRACT 

We study global inertial-modes with the purpose of unraveling the role they play in the tidal 
dissipation process of Jupiter. For spheres of uniformly rotating, neutrally buoyant fluid, we show 
that the partial differential equation governing inertial-modes can be separated into two ordinary 
differential equations when the density is constant, or when the density has a power-law dependence 
on radius. For more general density dependencies, we show that one can obtain an approximate 
solution to the inertial-modes that is accurate to the second order in wave- vector. Frequencies of 
inertial-modes are limited to to < 2f2 (51 is the rotation rate), with modes propagating closer to the 
rotation axis having higher frequencies. An inertial-mode propagates throughout much of the sphere 
with a relatively constant wavelength, and a wave amplitude that scales with density as 1/ \fp. It is 
reflected near the surface at a depth that depends on latitude, with the depth being much shallower 
near the special latitudes 9 = cos -1 ±cj/2fi. Around this region, this mode has the highest wave 
amplitude as well as the sharpest spatial gradient (the "singularity belt"), thereby incurring the 
strongest turbulent dissipation. Inertial-modes naturally cause small Eulerian density perturbations, 
so they are only weakly coupled to the tidal potential. In a companion paper, we attempt to apply 
these results to the problem of tidal dissipation in Jupiter. 

Subject headings: hydrodynamics — waves — planets and satellites: individual (Jupiter) — stars: 
oscillations — stars: rotation — convection 



1. INTRODUCTION 
1.1. Physical Motivation 

We study properties of global inertial-modes in neu- 
trally buoyant, rotating spheres. We briefly intro- 
duce our motivation for doing so, but refer readers to 
IQgilvie fc Lml <|2004fl . as well as[ffiiJ l|2004 hereafter Pa- 
per II) for a more thorough introduction. 

Jupiter's tidal dissipation factor (dimcnsionless Q 
valu e) has been estimated to be 10 5 < Q < 2 x 
10 6 JOoldreich fc SoteHlSlRTaTe fc Oreenherdrmall. 
based on the current resonant configuration of the 
Galilean satellites. Q is inversely proportional to the 
rate of tidal dissipation. So far, the physical explana- 
tion for this Q value has remained elusive, with the most 
reliable calculations yielding a tidal Q orders of mag- 
nitude greater than the above inferred value. Impor- 
tantly, known extra-solar pl anets also exhibits Q value 
of the same order as Jupiter ( Wu 2003) . This is inferred 
from the value of the semi-major axis (~ 0.07 AU) within 
which extra-solar planets are observed to possess largely 
circularized orbits. 

We call attention to two common characteristics shared 
by Jupiter and the close-in exo-planets: except for a 
thin radiative atmosphere, the interior of these bodies 
are convective; and they spin fast - typical spin periods 
are shorter or comparable to twice the tidal forcing pe- 
riod. This motivates us to search for an answer to the 
tidal Q problem that relates to rotation. In a sphere of 
neutrally buoyant fluid, rotation gives rise to a branch of 
eigenmodes: inertial-modes. These modes are restored 
by Coriolis force and have many interesting properties 
that make them good candidates for explaining the tidal 
Q value. 

1.2. Previous works 



Waves restored by Coriolis force have been studied 
extensively in the oceanographic and atmospheric sci- 
ences. Various names are associated with them, e.g., 
Rossby waves, pla netary waves, inertial waves, R-modes 
l)Oreenspanll 1 968|k In these contexts, the waves are typ- 
ically assumed to propagate inside a thin spherical shell. 

There have been a few studies of inertial-m odes 1 in 
the astrophysical context. iSchenk et alJ l)2002f) gives a 
fairly complete survey of the literature related to inertial- 
modes. We refer readers to that paper for a better un- 
derstanding of the nomenclature and past efforts. Here, 
we only mention a few related early works that have par- 
ti cular impact o n the current one. 

iBrvanl JT889) studied tidal forcing of oscillations in 
a rotating, uniform density spheroid (or ellipsoid). He 
performed a coordinate transformation (pioneered by 
Poincarc) under which the oscillation equation becomes 
separable. This immense ly facilitates ours a nd ot hers' 
study of inertial-modes. iLindblom & Ioseil l)1999l ) fol- 
l owed essentially the same app roach. 

Papaloiz ou fc Pringld l)1981|) studied inertial-modes in 
a fully convective star with an adiabatic index T = 5/3 
(equivalent to an = 1.5 polytrope). They realized that 
one could obtain the eigenfrequency spectrum fairly ac- 
curately without detailed knowledge of the eigenfunction. 
This is achieved using the variational principle and by ex- 
panding the eigenfunction in a well-chosen basis. They 
also pointed out that the fully convective case is special 
in that gravity-modes vanish so inertial-modes form a 
single sequence in frequency that depends only on rota- 
ti onal frequency. 

lLockitch &: Friedman! l)1999l LF from now on) studied 
inertial-modes in spheres with arbitrary polytropic den- 

1 Inertial-modes (also called generalized r-modes. hybrid modes) 
refer to rotationally restored modes in zero-buoyancy environment. 
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sity profile. They expressed the spatial structure of each 
inertial-mode as a sum of spherical harmonic functions 
and curls of the spherical harmonics. They presented 
some eigenfrequencies for, e.g., n = 1 polytrope. We will 
compare our results against theirs. 

Incrtial-modes have also been att acked numerically, via 
integ ration of the characteristics l|Dintrans Ez Rieutordl 
120001 lOeilvie fc Linl 120041) . th e finite difference method 
llSavoniie &: PaDaloizoulll99^K and the spectral method 
IQgilvie fc Linl i2004D . 

1.3. This Work 

There are two purposes to this paper. First, it lays 
down the foundation for Paper II where we discuss our at- 
tempt to solve the tidal Q problem. Second, it presents a 
new series of exact solutions to inertial-modes in spheres 
with power-law density profile, as well as approximate 
solutions for spheres with an arbitrary (but smoothly 
varying) density profile. 

Our approach centers on the ability to reduce the par- 
tial differential equation governing fluid motion in a ro- 
tating sphere to ordinary differential equations. This 
semi-analytical approach produces results that are both 
easily reproducible and have clear physical interpreta- 
tions. The mathematics are concentrated in <J21 We then 
expose properties of inertial-modes that are relevant for 
its interaction with the tidal perturbations (©. A dis- 
cussion section (<2J follows in which we consider the va- 
lidity of our assumptions. 

Readers interested in the tidal dissipation problem 
alone are referred to for a brief summary of the fea- 
tures of inertial-modes upon which we build our theory 
of tidal dissipation (Paper II). 

Readers interested in getting only a favor of how 
inertial-modes look like are referred to Figs. 0]&[5| The 
accompanying description f ij3.4l) as well as a comparison 
between inertial-modes and the more commonly known 
gravity- and pressure-modes ( H3.5| l may prove useful as 
well. 

2. INERTIAL-MODE EIGENFUNCTIONS 

This section documents our effort in obtaining semi- 
analytical eigenfunctions for inertial-modes. The rele- 
vant equation of motion is first introduced in M2.ll We 
then deal with increasingly more complicated and more 
realistic cases. This proceeds from simple uniform den- 
sity spheres ( ^2.2(1 . to power-law density spheres ( H2.3| l. 
and lastly to spheres with realistic planetary density pro- 
files (1231 

2.1. Equations of Motion 

Consider a planet spinning with a uniform angular ve- 
locity of H pointing in the z direction. In the rotating 
frame, the equations for momentum and mass conserva- 
tion read 

p£ + 2pftx£ = -Vf/+— p'-pV&, (1) 

P 



V •&*) = (), 



(2) 



where £ is the displacement vector, while p' , p' and $' are 
the Eulerian perturbations to pressure, density and grav- 
itational potential, respectively. We ignore rotational de- 
formation to the hydrostatic structure, as well as the 



centrifugal force associated with the perturbation. Both 
these terms are smaller by a factor ~ (ft/ft max ) 2 than 
the terms we keep, with f2 max being the break-up spin- 
rate of the planet. This factor is ~ 10% for Jupiter and 
much smaller for close-in extra-solar planets which likely 
have reached spin-synchronization with the orbit. 

We restrict ourselves to adiabatic perturbations, hence 
the Lagrangian pressure and density perturbations are 
related to each other by 



Sp Sp 
— — 1 1 — i 
P P 



(3) 



where the adiabatic index T\ = d\np/ din p\ s and is re- 



lated to the speed of sound by Ti 



'p/p. 



The interiors of giant planets are convectively unstable, 
with a low degree of super-adiabaticity, as is guaranteed 
by the fact that the convective velocity is fairly subsonic. 
This allows us to treat the fluid as neutrally buoyant. 
Setting the Brunt- Vaisala, frequency to zero, we obtain, 



dp p dP 
dr Tip dr 



(4) 



Given the following expression which relates the La- 
grangian and the Eulerian perturbations in a quantity 
X, 

5X = X' + £- VI, (5) 
equations © and Q combine to yield 



P F P 
— — 1 1 — 

P P 



(6) 



and the right-hand side of equation Q can be simplified 
into 



Vp' Vp p' 



P 



(7) 



(8) 



P P P VP. 
where we have introduced a new scalar 

lo z \ p J lo z \ p 

with lo being the mode frequency in the rotating frame. 
We adopt for all variables the following dependence 
on time (t) and on the azimuthal angle (4>): X cx 
exp [i(m<f> — u>t)]. Modes with denotation (m, to) and 
(—to, —lo) are physically the same mode, so we restrict 
ourselves to lo > 0, with m > representing a prograde 
mode, and m < a retrograde one. 

In solving for the inertial-mode eigenfunction, we 
adopt the Cowling approximation (negligible potential 
perturbations associated with the fluid movement) and 
ignore any external potential forcing (e.g., the tidal po- 
tential), so $' = 0. The former is justified in H4.ll and 
the latter is relevant as we are interested in free oscilla- 
tions. Equation |S| yields, 



/ ^ 2 P i 



(9) 



Following convention, we define the following two di- 
mensionless numbers: 



LO 

2ft' 



1 



(10) 
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Inertial- modes satisfy < fi < 1. Equations CQ-J2J) can 
now be recast as 



i + iq{e z x £) = Vip, 

u' 2 



(11) 
(12) 



where = —dr/d\np is the density scale height and 
H = Cg/g (eq. 0]), with 3 being the local gravitational 
acceleration. 

Operating on equation 1)11(1 with e z ■ and e z x , respec- 
tively, and replacing e z ■ £ and e z x £ using the resultant 
equations, we obtain the following relationship between 
£ and ip 

£ = —L—(l-iqe z x -q 2 e z e z -) Vip. (13) 
1 — q z 

We substitute this equation into equation l|12|) to elimi- 
nate £ and to acquire the following key equation for tp, 



V 2 V - q 2 

L (f*t 
H \~dr~ 



d 2 ip 
dz 2 



q cost 



^dip mq 
dz 



u> 



^ -(l-g^)^.(14) 



Here 9 is the zenith angle and cos# = z/r with z be- 
ing the height along the rotational axis. Let w be 
the cylindrical radius. The partial derivatives here are 
to be understood as d/dr = d/dr\ e , dj 89 — d/88\ r , 
d/dz = d/8z\ m , and d/dm = d/dw\ z . In general, the 
above partial differential equation is not separable in any 
coordinates and only fully numerical solutions could be 
sought. 

The displacement vector in spherical coordinates is re- 
lated to ip as lfLl|l . 

1 / dip mq 2 dip\ 

^W^-T^ 5 C0S ^J' (15) 

C~(-^-^ + , 2 sm^), (16) 
1 — q^ \r 09 rm dz J 



If im dip 

£0 = -, 2 — tp - iq^- 

1 — q^ \ m dm 



(17) 



In the following discussions (except where noted), we 
scale all lengths by the radius of the planet (R). 

Equation l|14|) describes the pulsation for both inertial- 
modes and pressure-modes in a uniformly rotating, 
neutrally-buoyant fluid - gravity-modes have zero fre- 
quencies in such a medium. While the compressional 
term ((1 — q 2 )uj 2 / c 2 ip) represents the main restoring force 
for the pressure-modes, it is negligible (in comparison to 
the Coriolis terms) for inertial-modes as the latter are 
much lower in frequency. We adhere to this simplifica- 
tion throughout our analysis, and further justify it in 

S3 

We now introduce an important coordinate tra nsform 
which facilitates much of our analysis. We follow IBrvanI 
(1889) in adopting a set of ellipsoidal coordinates (xi, X2) 
which depend on the value of fi. More details on these 
coordinates are presented in Appendix Fig. 2] de- 
picts the topographic curves of (x\,X2) in a meridional 
plane, as well as the curves of constant r in the X\ — x-i 
plane, when fi — 0.75. Generally, the x\ coordinate is 
parallel to the m axis, while the X2 coordinate is parallel 



to the z axis over much of the sphere. On the spherical 
surface, either x\ or X2 (or both at the special latitude 
9 = cos -1 ±/i) equals fi; x\ = 1 at the rotation axis and 
x% = at the equator. 

The ellipsoidal coordinates are a natural choice for 
studying inertial-modes in a rotating sphere Motion re- 
stored by the Coriolis force (inertial-waves) would have 
liked to follow the cylindrical coordinates, were it not for 
the spherical topology of the object within which they 
dwell. The ellipsoidal coordinates, being a hybrid be- 
tween the cylindrical and the spherical coordinates, are 
therefore particularly suitable for this purpose. 

With this set of coordinates, the left-hand side of equa- 
tion H14fl is transformed into, 



i 2 ^ ) V; 



I 2 

dz 2 , 

1 d ( d \ Id 2 2 . d 2 

m dm \ dm J m l dcp z dz z 



1-A* 2 



d 2 d 



d 2 d 

(1- X 2)^2-2X2 — 



' dx 2 



1 



ip. 



(18) 



l 2 u^2 — ->-2 . 

This, as we shall see shortly, allows the separation of 
variables under certain circumstances. 

2.2. Uniform Density Sphere 

2.2.1. Formal Solution 

In a uniform density sphere, the scale height H = 00 
and the right-hand side of equation ((14(1 vanishes. Its 
left-hand side (eq [18] ) can benefit from the following 
decomposition, 



ip = ipi{xi) ip 2 (x 2 ) 



with ipi satisfying 

ViiPi + K 2 tp t = 0. 
Here the differential operator T>i is 



Vi = 



d 

dxi 



(1 



I — 

dxi 



(19) 



(20) 



(21) 



and if is a constant introduced when we separate vari- 
ables. 

I ■ result was first obtained bv IBrvanI l|1889|) and its 
solutions are called 'Bryan's modes'. In fact, the solu- 
tions to ip\ and ip2 are the associated Legendre poly- 
nomials. Requiring ip to be finite at the rotation axis 
(x\ = 1), we find ipi and ipi to be the same spherica l 
harmonic of the first kind ((Abramowitz fc Stegunl ll972'l: 
ipi = i>2 = PPix) with £ being an integer, K 2 = l{t + l), 
and the variable x taken over the ranges x\ G [fi, 1], and 
x 2 e [— fi, fi], respectively. We explicitly require that 
ipi (x 1 = fi) = ip2{x 2 = fi) so the eigenfunction needs 
only one normalization constant. 

The following boundary conditions apply. First, at the 
equator {x% — 0), even-parity modes 2 satisfy 



dip2 

d,X2 



0. 



(22) 



x 2 =0 



2 It is straightforward to show that the displacement vector £ 
has the same equatorial symmetry as -0. 
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while odd-parity modes satisfy 



2\x 2 =0 



0. 



(23) 



Properties of the Legendre polynomials 

ijAbramowitz fc Stegunl 11972(1 require that (£ + to) 
to be an even integer in the former case, and odd in the 
latter. Second, ipi is finite at the polar axis (x\ = 1). 
The numerical equivalent to this statement is best 
realized by introducing a variable gi which is related to 
ipi as 

^( Xi ) = (1 - x?)H/2 ff .( x .). 



This variable satisfies (eq. |2()(l 
d 2 g 



(1-4) 



dx] 



2x41ml 



axj 



0, 



(24) 



(25) 



where A 2 = K 2 - \m\(\m\ + 1) = £{£+l) - \m\(\m\ + 1). 
Regularity of the eigenfunction at x\ — 1 then translates 
into a boundary condition 



dgi 

dx\ 



X 2 



2(\m\ 



1 



-01 • 



(26) 



a:i — 1 

Our solutions show that near the rotation axis, g\ ap- 
proaches a constant, while ipi approaches zero (if |m| > 

In this problem, the eigenfunction can be solved in- 
dependently of the eigenvalue /i. So to determine /i, we 
need one more boundary condition. We enforce the phys- 
ical condition that there is vacuum outside the planetary 
surface (r = 1), and therefore pressure perturbation at 
the surface h as to be zero. Wr itten in a convenient form 
of dp/p = ||U nno et aLll!989(l . this corresponds to (us- 
ing eq. |12j and ignoring the compressional term), 



dp 
P 



Sp^p 

p p 



-ri(v-o 



p 



-gCr = 0. 



(27) 



We relate £ r to ip using equation (|15(l , as well as relations 
presented in Appendix 1X1 



1 



(1 



p 2 ){x\-x 2 ){l 
dip [i 2 — x\ 
dxi 1 



Xl 



x 2 



X2 



x\){l-xl) 

dtp [A 2 — 2 



8X2 1 



+TO/i 



(l-x 2 )(l-x 2 ) 



So requiring £ r = at the surface {x\ 
is equivalent to requiring that 



(i or \x 2 \ 



(28) 
= *0 



#1 
dx\ 

dip2 
dx 2 



x 1= fi ■ 



-SIGN[ar 2 ] 



\x 2 \=fi 



1 



1p2 



\x 2 \=ji 



(29) 



Since ipi and ip 2 have the same functional form (P £ m (;r)), 
these two equations are in fact the same thing. In actual 
numerical procedure, we solve for gi as opposed to ipi. 
Equation i|29[l is modified for gi as, 



dgi 
dxi 



— (to — |m|/i) 
1 -/i 2 



9i\x 1=t i- 



(30) 



3 For each (£, to) pair, there can there- 
roots satisfying this boundary condi- 



Using equation 1)24(1 . we find that gi is a polynomial of 
order (£ — |m| 
fore be (£ — \m\ 
tion, with half of them having /i > 0. Another way of 
phrasing this is that, for each eigenfunction {ip{x\, x 2 ) — 

-Pj" l '(xi) Pj m '(x2)), there are (£ — |m|) eigenfrequencies. 
Roughly half of these are prograde modes (to > 0) and 
the other half retrograde modes (m < 0). 

Table ^ presents some example eigenfrequencies for 
inertial-modes in a uniform density sphere. Our results 
agr ee with those presented in Tables 3-4 of LF and Table 
1 of lLindblom fc Ipserl <|1999|1 . 

2.2.2. The Dispersion Relation 

In the previous section, we have introduced various 
eigenvalues (K, A, £ and to), as well as the eigenfrequency 
li. What is the geometrical meaning of these eigenvalues 
and what is the dispersion relation for inertial-modes? 

Here, we present a derivation that answers these ques- 
tions. We first convert the independent variable from xi 
to @i = cos -1 Xi (to be differentiated with the spherical 
angle 9) so that equation f2*5|) becomes, 



d 2 g l , , , , . ^cosO, dgi 2 
—2 + (2|to| + 1)377^— + A g % = 0. 



d&j 



sin Qi dQi 



(31) 



Adopting a WKB approach, we express gi cx 
exp(z J k^dQi) with the wave vector fee = k^+ikj where 
the real part kn ~ 0(A) ^> 1 and the imaginary part 
k] ~ 0(1). Substituting this into equation (|31|) . and 
equating terms of comparable magnitudes, we find 

kn~ A, 

/, l\ cosOi /, l\ din sin 9; .„„, 

fc K H+ 2Wn H+ 2j ^e— (32) 

Together with equation l|24(l , this gives rise to the follow- 
ing approximate solution for ipi in the WKB regime, 



iPi = pr( Xl ) cx 



sin Oil 1 / 



YH cos(A6, + a) 



(33) 



So ipi is an oscillating function with a roughly constant 
envelope. Moreover, a = — Xn/2 + Tr for even-parity 
modes (eq. [22); while a = — \n/2 + n/2 for odd-parity 
modes (eq. [231 ) - Nodes of ipi = are roughly evenly 
spaced in 9 e [0,7r/2] space, and the value of A cor- 
responds with twice the total number of nodes (more 
below). 

To obtain the dispersion relation, we insert the above 
approximate solution into the boundary condition at the 
surface (eq. [21]) where 9q = cos -1 ±/.i. This yields, 



Asin(A9o + a) + 



For A> 1, this can be approximated by 



cos(A9 + a) w 0. (34) 



sin(A9o + a) m A9q + a 



±TO 



HIT : 



(35) 



where n is an integer. We retain only the positive sign 
as to can be positive or negative. For even-parity modes, 

3 One can express gi as d) m \ P®(x)/dx\ m \ , and P®(x) is a poly- 
nomial of order I. 
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the dispersion relation for the eigenfrequency runs as 



_ / nir m 

p = cos Bo = cos — — 



A V /T^2 A 2 2 A^ 
(1 + n)ir m \ 



while for odd-parity, it is 



(36) 



= COS — ■ 



A v/T^A 2 



7T 

2 



(i +n)7r 



A 



v/W^A 2 



(37) 



If we define the number of nodes in the 6j S [0, 9 ] 
range to be n 1; and the number in the 2 € [©o,7r/2] 
range to be ri2, they satisfy ni7r w A8o, n2^ ~ \(ir/2 — 
0o). So A w 2(ni + 712), as we have mentioned ear- 
lier. Moreover, equation (|35|) yields n f=a n2 — 1 + 
m/ y/l — p 2 \ir for even modes, and n w 712 — 1/2 + 
m/y/l — fjPXir for odd modes. Substituting these esti- 
mates into equations lj3U|> & (|3TJl . we find the following 
simplified dispersion relation, 



to 

m 



n 2 TT 



2{m + n 2 ) 



2{n\ + 122) 



( 38 ) 

This dispersion relation has also been numerically con- 
firmed. So for a given set of (A,m), the eigenfrequency 
p rises as the partition of the x\jxi space moves from 
p = (when m ~ A/2, n2 = 0) to p = 1 (when m = 0, 
n2 = A/2), and there are ~ A/ 2 discreet eigenfrequen- 
cies. Recall that A = y/HJ+T) - \m\(\m\ + 1) w i- \m\. 
So for each (A, \m\) pair, we find ~ (I — \m\)/2 prograde 
eigenmodes and ~ (£— \m\)/2 retrograde ones, consistent 
with the results found in H2.2.1I 

In H3.ll we will discuss further the geometrical meaning 
of the eigenvalues in spherical coordinates. 

2.3. Power-Law Density Sphere 
2.3.1. The equation is Separable! 

We are inspired by the planetary density profile to in- 
vestigate inertial-modes in a sphere where density de- 
pends on radius as a power-law. To our pleasant sur- 
prise, the equation turns out be separable in this case. 
This forms the basis for much of our work when dealing 
with real planets. 

We adopt the following power-law density profile (see 
eq. HH) 

p ex (1 - r 2 Y j = po [{x\ - p 2 ){p 2 - xl)] P . (39) 

This is not a usual choice and is different from the more 
typically discussed polytropes. However, near the sur- 
face, our power-law density profile behaves like a normal 
polytrope. Let S = 1 — r be the depth into the planet. 
For (J <C 1, we find 



5 = 1 



1 



0! - P 2 )(P 2 



A) 



(1 - P 2 )p 2 



IM. 

2 



p 2 ){p 2 -xl) 



(1 



p 2 )p 2 



(40) 



So near the surface, equation l(3^|) yields p cx S 13 , as is 
the case for a polytrope model with the polytropic index 
n = j3. The advantage of choosing such a profile will 
become clear soon. 

We obtain the following expression for the density 
scale-height H: 

dhip fzd^wd 
dr \r dz r dm 

(1 - p 2 )p 2 (3r 



fi- 



llip 



( x i — p 2 )(p 2 — A) ' 

where equ ations (IA2D & l|A6(l are used. Together with 
equations (|A4() - (|A6() . this allows equation l|14|) to be re- 
cast into (again ignoring the comprcssional term), 



T>2 + 



2(3 Xl (l~x 2 ) d 
x\ — p 2 dx\ 

2(3x 2 {l-xl) d 
x\ — p 2 8x2 



2p(3m 



x\ - p* 
2p[3m 



iP = 0, (42) 



where the operator T>i is defined in equation l|21|) . So 
when density follows the power-law profile (eq. [3D]), 
the equation for the inertial-modes is again separable: 
t/j(xi,x 2 ) = i>i(xi)ip2(x2) with fa satisfying 



2(3 Xl {\ 



x- - p' 

*2„ 



dxi 



2p(3m 
x 2 - p 2 



£nl>i + K*i/>i = 0, (43) 

where we introduce a new operator £i. Note that this 
equation is exact (except for omitting the compressional 
term) for the power-law density sphere. 

Numerically, it is more accurate to solve for g, = 
ipi /(I — x 2 )'" 1 '/ 2 . This function satisfies 



(1-*?) 



A 



d 2 9i 
dx\ 
2(3\n 



2x l (\m\ + 1) 
x 2 2(3mp 



dgi 2(3x t {l - x 2 ) dgt 
dxi x 2 — 1 



dxi 



x 2 - p 2 



P z 



9i = 0, 



(44) 



where A 2 = K 2 — lm|(|m| + 1). The boundary conditions 
presented in ij2.2.1l applv here as well, with the exception 
of equation (|26|l which should be modified into 

dgr = A 2 + 2(3^ 

dx lxi=1 2(|m| + l) 91 ' [b) 

So unlike the uniform density case, both the inertial- 
mode equation and the boundary conditions now depend 
explicitly on the value of p as well as the sign of m. 
This breaks the eigenfunction degeneracy existing in the 
uniform density case: each eigenmode now has a unique 
eigenfunction. 

In Fig. |2 we plot the resulting eigenfunction for a 
m = —2 mode in a (3 = 1 model. We also contrast the 
density profile of our (3=1 model with the conventional 
n = 1 polytrope model (p cx p 2 ). Details of how we 
construct the (3 models are presented in Appendix iBl 

Table El lists some sample eigenfrequencies for \m\ = 2 
modes in various (3 models. The number of eigenmodes 
remains conserved when (3 varies, with a close one-to- 
one correspondence between modes in different density 
profiles. This makes mode typing trivial. We also com- 
pare our results against those of LF obtained for a n = 1 
polytrope. 
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2.3.2. WKB Envelope 

While the eigenfunction ip inside a uniform density 
sphere has a roughly constant envelope (eq. \6'6\ ) . ip in 
the power-law density case behaves differently. Following 
the same operation as in W2.2.2\ one finds 



Introducing a variable t, 

t=(x\- p 2 )(p 2 -x\) 
one finds, 

dlnX 



p 2 (l-p 2 )(l-r 2 ), (50) 




2 2f3(mp — \m\ cos 2 9,) 

A ~\~ i r, ~ ~ ~ A 



dr 



-2rp\\-p 2 ) 



dlnX 
dt 



(51) 



1 

m\ + - 



j cos 2 8i — /i 2 | 
cos Qi /3cos9isin9i 



Now repeat the same calculations that lead to equation 
<{4*2l , we obtain 



and 



sin Qi 
lnsin |m|+ 5 9, 

1 



cos Of — p 2 



In 



9 ? , 



,(46) 



(Si - £ 2 )V> 
, (l 



dlnX 



dint 



(1 - xl)xi d 



(x\ — p 2 ) dxi 



x\)x 2 d 



mp(x 2 — a; 2 ,)" 1 



(p 2 - x\) dx 2 



*/> = 0, (52) 



I sin 9 ? - 1 2 I cos 2 9, 



■cos(A9i+a). (47) 



So, the eigenfunction ip has an amplitude that rises to- 
ward the surface (when \xi\ — » p). In fact, its envelope 
scales with density as 



where the operator £i is defined in equation (|43|l . Here, 
since dtp/dxi ~ Xip, we can not ignore terms in the 
square parenthesis if we want to be accurate to 0(1/ A 2 ), 
the final goal of our procedure. 

Instead, we experiment with the following decomposi- 
tion for ip, 



1p = i>l{x 1 )lp 2 (x2) OC 



1 



1 



ip = VXip (xi,x 2 ) = 



where 



dVx 

dxi 



(£ x - £2)^0 + (d x 
where the coefficient d* is 



[(x\ - p 2 ) (p 2 - x 2 )]P/ 2 VP' 

(48) 

Fig. shows one example of such a WKB envelope. 

We find that the dispersion relation is also modified 
slightly from the case of a uniform density sphere. But 
qualitative features of equation l|38f) remain. 

2.4. Approximate Solution for Realistic Density Profiles £iip = \^X£itpo J r'^ a i 

The density profile inside a planet typically traces out 
two different polytropes: near the surface, the gas can 
be approximated as an ideal gas composed of diatomic 
molecules, with a mean degree of freedom 5, so Ti = 
dlnP/dhip\ s — 7/5, and the correct polytrope number 
is n ~ 2; in the interior of the planet, Coulomb pressure 
and electron degeneracy modify the equation of state and 
raise T\ to ~ 2, while reducing the polytrope number 
to n ~ 1. This motivates us to model Jupiter's density 
profile using two power-laws of the form in equation : 
(3 = 1 in the interior and /3 = 1.8 near the surface, with 
the transition occurring at a radius r ~ 0.98 (see Paper 
II). This reproduce s profiles in realistic Jupiter models 
JGiiillot et alJ l2004). The presence of a core or a phase 
transition complicates this picture. We discuss them in 
more detail in Paper II. 

In general, equation 114|) is not separable under these 
density profiles. However, we discover an approximate 
solution which is accurate to the second order in wave- 
numbers (C(l/A 2 ) with A 3> 1), for the case when 
the density profile is power-law near the surface and 
smoothly varying in the interior. This is largely inspired 
by the WKB result in WH\ 

We first introduce a fiducial density jO sur f which satisfies 
equation (|39|l with a power-law index taken to be that 
for the true density p near the surface. We also define 
X = psuil/p so X fa 1 near the surface and deviates 
from unity toward the center. Inverse of the density scale 
height can be written as 



Psurf 



tpo(xi,x 2 ) 



(53) 



This is inspired by the WKB envelope presented in equa- 
tion J2HJ). 

Formally expressing £iip = aid 2 ip / dxf + bidip/dxi + 
Ciip, where the meanings of a,, 6, and Cj are clear from 
equation (|43|l . we find 



dyX chj)g 
dxi dxi 



d 2 Vx , dVx 



dx 2 



dlnX 
dint ' 



dxi 



(54) 
(55) 



A straight-forward but lengthy derivation shows that 
equation (|52|l can be recast into an equation for ipQ, 



d 2 )ip = 0, 



(56) 




Xi(l-x$)dVX 



(1-x 2 )^ 

[ J dxi 

2mp 



2(3-2 



dlnX 
dint 



dlnX 



dx-. 



p 2 dint 



(57) 



1_ dlnp _ dlnX 2(1- p 2 )p 2 f3r 
H dr dr (x\ — p 2 )(p 2 — x?,) ' 



(49) 



Here refers to the power-law index for p S urf, or that 
for p near the surface. The benefit of the transforma- 
tion introduced in equation i]")-!!) is that di contains no 
derivatives on ip. 

Near the surface, X = const, d t — as is expected 
and equation (|56|l is separable. But in deeper region of 
the planet, X is a complex function of X\ and x 2 and 
equation (|56J) is not separable. 

However, if X is a smoothly varying function with a 
scale length being the radius of the planet 4 , one can show 
that, in the WKB region, diipo ~ ipo and is ~ 1/A 2 
smaller than the fyipo term. So if we ignore the di terms, 

4 In the case of Jupiter, X = 1 near the surface and gradually 
rises to larger values toward the center. 



Tides in Jupiter 



7 



we only introduce errors of order 0(A~ 2 ) <C 1. It is worth 
pointing out that the non-WKB region occurs near the 
surface where our solution is exact. 

In conclusion, we can adopt the following approximate 
solution for tf> 



Psurf 



4>1 (xi)lb 2 (x 2 ), 



(58) 



with ipi(xi) satisfying Eiipi + K ipi = 0. This solution 
is exact near the surface and is accurate to 0(1/ X 2 ) in 
the WKB region. The latter attribute indicates that the 
approximate solution describes accurately both the en- 
velope and the phase of the actual inertial-mode eigen- 
function. 

3. PROPERTIES OF INERTIAL-MODES 

Having discussed methods to obtain inertial-mode 
eigenfunctions in various density profiles, here we focus 
on studying general properties of inertial-modes. This 
include the WKB properties, the normalization relation- 
ship and density of modes in the frequency spectrum. 
Moreover, we attempt to give readers a graphical im- 
pression of how inertial-modes look like both inside and 
outside the planet. Lastly, we compare inertial-modes 
against well known gravity- and pressure-modes to gain 
intuition for this branch of eigenmodes. 

3.1. WKB Properties 

We first derive the WKB dispersion relation for 
inertial-modes. We then determine the confine of the 
WKB region, and end with a general derivation for the 
WKB envelope of inertial-mode amplitude. 

In the WKB region, let V m ik and d/dz « ik z , equa- 
tion (jTl|) yields, 

k 2 - q 2 k 2 « 0, (59) 

or /i pa k z /k < 1. This is to be compared with result 
from a more careful derivation (' H2.2.2|l which shows /i w 
sin(n 2 ir/X). Since the x 2 axis is largely along the z axis 
(Fig. [TJ, k z ~ n 2 /R. So we have k ~ X/R ~ 2(ni + 
n 2 )/R, with R being the planet radius (normalized to be 
one). Such a dispersion relation implies that the mode 
frequency (lo = 2Qfi) does not depend on the number of 
wiggles in a mode, but rather on the direction of wave 
propagation. Modes that propagate close to the rotation 
axis have higher frequencies (lo ~ 2ft) than those that 
propagate close to the equator (lo ~ 0). 

Such a dispersion relation can be understood by the 
following physical argument. We relate £ to tp using 
equation ifblf . 



, fc - q 2 e z k 2 
1-q 2 



(60) 



So the fluid velocity (v — d£/dt) is perpendicular to the 
phase velocity (u p h = Lo/k 2 k), or v ■ v p h ~ 0. Inertial- 
waves are largely transverse waves. A mode with its 
phase velocity along the equator will show fluid motion 
in the z direction. As a result, it experiences little Cori- 
olis force, and has a frequency that is close to zero. In 
contrast, a mode with k along the z direction will ex- 
perience the strongest restoring force and will have the 
highest frequency. 



While the phase velocity of the inertial-wave ■Uph : 
lo /k 2 k, the group velocity of the inertial-wave runs as 



duo 
dkz 



dkh 



e h = 2VL 



k 2 
fc 3 



-e h 



:(k - q 2 e z k z ), 



(61) 



where k = k z e z +kheh- It is easy to show that v p h -v g = 
0. 

Now we search for the boundary that separates the 
WKB cavity from the evanescent cavity for inertial- 
modes. Consider propagation in either one of the (x\ , x 2 ) 
directions. Equation l|46|) shows that the real part of 
the wave-vector remains fairly constant over the whole 
planet, and is little affected by the density profile, while 
the imaginary part of the wave-vector rises toward the 
surface as (3/\xi — p\. The latter becomes comparable 
to the former when \xt — fJ,\ < (3/X. Recalling that 
\xi\ = [i occurs at the surface, one sees that an upward 
propagating wave is reflected near the surface at a depth 
S ~ /3/X at most latitudes, except near a special latitude 
| cos#| ~ fj, when the wave penetrates much higher into 
the envelope, to a depth of S ~ (3 2 /X 2 . We call this spe- 
cial latitude where xi ~ \x 2 \ ~ /x the "singularity belt" 
and it is an important region for mode dissipation. We 
will return to this concept in ^13 .41 

Lastly, we turn to study the WKB envelope of an 
inertial-mode. Multiplying equation (JJJ by £, and sim- 
plifying the resulting expression using equations 10} & 
© , we arrive at the following equation of energy conser- 
vation, 



9t I 2 ? 



P 



2pc] 



V • (jp'k) = o. 



(62) 



Terms in the first set of parenthesis are readily identi- 
fied as the energy density (including both kinetic energy 
density and compressional energy density), while that in 
the second set is the energy flux carried by the wave. 
Both quantities look identical to those for non-rotating 
objects. This is expected since the Coriolis force is an 
inertial force and does not do work or contribute to en- 
ergy. 

The relative importance between the two energy den- 
sity terms is, 

^PM? c 2 s \e cjA^ ^ / c 2 X 2 \ ( GM/RT 
j 2 R 2 ~ \GM/RJ \ lo 2 



P' /PCs 



LJ 2 1p 2 



(63) 

Since we are interested in A > 1 modes, and since 
for planets rotating well below the break-up speed, 
inertial-mode frequency (lo < 2Q) is smaller than the 

fundamental frequency of the planet ( yj GM /R 3 ), the 
above ratio is much greater than unity over much of 
the planet. Exceptions occur near the surface, where 
S < 1/X 2 (lo 2 /GM/R 3 ) < 1/A 2 , well outside the WKB 
cavity. So the energy density of an inertial-mode is domi- 
nated by its kinetic part. For a standing wave, the energy 
density (and in this case, the kinetic energy density) re- 
mains constant in time. So different velocity components 
must be out-of-phase with each other. For instance, near 
the surface, ~ \^\ » |£ r |, and £g and ^ are 90deg 
apart in phase. 
Take a standing wave of the form 

Z{r,t)=£{r)e- iut + ?{r)e iut , (64) 
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equation (|60|l gives rise to the following energy density 
(the time-independent part) 



k z -2q 2 k 2 z +q 4 k 2 



k z 



q 2 -l 



(l-q 2 ) 2 



(65) 



while the energy flux 

F=p'£ = -2ulm\p'(r)t(r)\ = -2w 3 p Im[ip £*] 
3 k - q 2 e z k 



q 2 - 1 



-p\i/>{r)\* = v g H, 



(66) 



where v g is the local group velocity as derived in equation 
<|61l) . Inside the WKB propagating cavity, the energy flux 
is constant (no reflection). This, coupled with the fact 
that k is largely constant inside the planet, yields the 
following WKB envelope 



Mr)\ 



(67) 



This agrees with results derived from more detailed con- 
siderations (eqs. |55 ] . H Q fc 

3.2. Mode Normalization 

We derive a scaling for the total energy in a mode. 
Applying results from ^3. H and the Jacobian defined in 
equation (|A3|) . we express total energy in a mode as 



E-- 



d A rU 



UJ 

T~q^ 

u) 2 p 



d 3 rp\ti(r)f 



d 3 rk 2 p\iP(r)f 



, | ^ 2 ^ 2 / {x\- x\)dxidx 2 d(j)k 2 p\ip\ 2 . (68) 

Most of the mode energy density lies inside the WKB 
region. Within this region, each nodal patch (AxiAo^) 
contributes a comparable amount to the total energy. 
Since the total number of patches is ~ m x n 2 , we obtain 



E' 



j^yn in2 (k 2 p\m\:z a x1 a X2 



X 2 ~0 



us 2 p 



{l-p?f 



nin 2 (p|^| 2 )|^~° 



(69) 



where the term p\ip\ 2 is supposed to be evaluated at 
the last crest away from x\ = 1 and x 2 = (center 
of the sphere). Adopting the arbitrary normalization of 
ip2(x2 = 0) = 1 and assuming that ipx(xi w 1) does 
not depend on n%, E oc n\n 2 if evaluated at a fixed fre- 
quency. Our numerical results, however, show (Fig. |3J) 
that E oc n 2 / 65 in a f3 = 1 model when only n% = n 2 
modes are considered (or p ~ 0.70). 

The small difference results from the fact that as ni 
increases, the amplitude of the last crest (ipi(xi « 1)) 
increases slightly. We adopt the numerical scaling E oc 



3.3. Density of Inertial- Modes 

For our tidal problem, it is useful to study how dense 
inertial-modes are within a given frequency interval. 
More specifically, we ask, for a given frequency po, now 
close is the closest resonance (smallest |/x — Moj) one can 
find among modes satisfying wave- number A < A max . 

Since frequencies of inertial-modes do not depend on 
the mode order, but on the direction of mode propaga- 
tion, modes of very different wave- numbers can coexist in 
the same frequency range. Fig. shows how p depends 
on A for even-parity, \m\ —2 modes in a uniform den- 
sity sphere. The mean frequency spacing between modes 
of the same A (or m + n 2 for more general models) is 
dp/dn 2 s» yl — p 2 ^/ A. This value decreases as p ap- 
proaches 1 or as A rises. The typical distance away from 
a resonance (\p — po\) is half this value. When we include 
all modes with A < A max , the closest resonance likely has 
\p - fi \ w y/l- p 2 n/\ 2 nax - tt/A: 



2 

max " 



,2.65 



for our later studies. 



3.4. How do Inertial- Modes Look Like? 

In this section, we provide a graphical impression for 
how inertial-modes look like both inside and on the sur- 
face of planets. This may be helpful for readers as 
inertial-modes are drastically different from modes in a 
non-rotating sphere, which can be described by a prod- 
uct of a radial function and a single spherical harmonic 
function. 

In the interior of a planet, an inertial- mode pro- 
duces alternating regions of compression and expansion, 
much like gravity- or pressure-modes do, except that 
for inertial-modes, these regions are lined-up along the 
(xi,x 2 ) coordinates. Fig. 0] depicts how the Eulerian 
density perturbation (//) and the perturbation velocity 
in the rotating frame look like in a meridional plane, for 
an example inertial-mode (ni = 5,n 2 = 3 and m = —2). 
There are three noteworthy features. The first is that 
the largest perturbations, as well as the steepest spatial 
gradients in these quantities, are to be found near the 
surface, especially near the angle | cos 9\ « p. The second 
feature is that velocity near the surface is purely horizon- 
tal. The radial component vanishes as is required by the 
boundary condition (eq. |27|V The third feature is that 
the velocity patterns inside the planet take the form of 
vortex rolls. Inertial-modes produce largely incompress- 
ible, and largely rotational motion (|V x £| 3> |V • £|). 

Fig. |U shows a surface view of the density perturba- 
tion for the same mode. For this retro-grade mode, the 
pattern rotates retrogradely on the planet surface. First 
notice the number of nodal patches on the surface. For 
a p- or g-mode in a non-rotating star, the number of ra- 
dial nodes does not show up in the surface pattern. For 
inertial-modes, however, the values of ni, n 2 and m (and 
even p) are all clearly embedded in the surface pattern. 
The surface pattern tells all. Second notice the presence 
of a belt (in both hemispheres) near | cos#| = p where 
the mode exhibits both the largest perturbation as well 
as the largest gradient of perturbation (see also Fig. 0J. 
This we call the "singularity belt" and it is a feature 
unique to inertial-modes. It is located at where both 
X\ — p and p — \x 2 \ < 1/A, corresponding to a region 
with a depth 5 ~ 1/A 2 and an angular extent 88 ~ 1/A. 
This region will turn out to be very important for tidal 
dissipation (Paper II). 
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For comparison, we present a meridional look for the 
m = — 2 R-mode in Fig. El R- modes are a special branch 
of inertial-modes, and we discuss them in H4.2I 

3.5. Comparison with Gravity- and Pressure-modes 

Modes restored by pressure or buoyancy (p- or g- 
modes) are familiar to astronomers. To help under- 
standing inertial-modes, we capitalize on this familiarity 
by discussing the differences between these modes and 
inertial-modes. 

Inertial-modes are more analogous to g-modes than p- 
modcs. Frequencies of inertial-modes are higher if their 
direction of propagation is more parallel to the rotation 
axis (ui/20, ~ k z /k, z being the rotation axis). Their fre- 
quencies are independent of the magnitudes of the wave- 
vector and are constrained to lo < 20. Similarly, higher 
frequency gravity-modes propagate more parallel to the 
potential surface (lu/N ~ kh/k, h being the horizontal 
direction) satisfying oj < N. Both inertial-waves and 
g- waves are transverse waves, i.e., their group velocities 
(direction of energy propagation) are perpendicular to 
their phase velocities. P-modes, in comparison, experi- 
ence stronger restoring force and hence higher frequen- 
cies if their wavelengths are shorter. Phase velocities and 
group velocities of p-modes are identical. 

Due to their low-frequency nature, both g-modes and 
inertial-modes cause little compression in their propagat- 
ing cavities. As a result, the energy of g-modes is domi- 
nated by the (gravitational) potential and kinetic terms, 
which are alternately important for half a cycle, each 
contributing in average half to the total energy. Inertial- 
modes live in neutrally stratified medium, their energy 
is dominated by kinetic energy alone. As a result, while 
g-modes (and p-modes) suffer dissipation from radiative 
diffusion and viscosity, inertial-modes are only sensitive 
to viscosity. 

As p- or g-waves propagate toward the surface, their 
wavelengths typically shorten, resulting in the largest 
dissipation. An inertial-mode, however, propagates in 
its WKB cavity with a roughly constant wavelength, ex- 
cept near the singularity belt (r ~ R and 9 ~ cos -1 ±/i) 
where its wavelength shrinks drastically. This is where 
we expect the largest dissipation to occur. 

Lastly, each p- or g-mode in a non-rotating star has an 
angular dependence that is described by a single spheri- 
cal harmonic function, Pf h {9, 4>). In contrast, the angular 
dependence of each inertial-mode is composed of a series 
of spherical harmonic functions. This implies that while 
for p- and g-modes, only the I = 2, \m\ = 2 branch can 
be excited by a potential force of the form P| (such as 
the lowest order tidal force), for inertial-modes, every 
\m\ = 2 even-parity mode can be excited. In this sense, 
the frequency spectrum of inertial-modes is dense, and 
the probability of finding a good frequency match (forc- 
ing frequency ~ mode frequency) is much improved over 
the non-rotating case ( H3.3p . 

4. FURTHER DISCUSSION 

4.1. Justifying Assumptions 

In our effort to obtain a semi-analytical solution for the 
inertial-modes (§2J, we have made a string of simplifying 
assumptions. We justify them here. 

In the analysis throughout this paper, we have ignored 
the compressional term (( 1 — ) w 2 /c 2 ^ in eq. ^3] ) . This 



is also called the 'anelastic approximation' (ignoring p'), 
and is often adopted when studying sub-sonic flows in 
stratified medium. This assumption is justified by equa- 
tion 1|63|) and discussions around it, where we show that 
the inertia term dominates over the compressional term 
in most of the planet except well above the WKB prop- 
agating region. Hence the compressional term does not 
significantly affect either the frequency or the structure 
of an inertial-mode. 

Inertial-modes are excited by the tidal potential 
through its density perturbation (p' , see Paper II). Could 
we be removing tidal forcing of inertial-modes by ig- 
noring p'l Fortunately, no. Equation © states that 
p' = pw 2 /c 2 ip. Since ip is non-zero, tidal forcing is zero 
only when c 2 — > oo. Taking the anelastic approximation 
does not preclude tidal coupling. It does imply, how- 
ever, that tidal forcing of inertial-modes is expected to 
be weaker than tidal forcing of the fundamental mode, 
by the same ratio that relates the inertia term to the 
compressional term. 

We have also ignored the potential perturbation caused 
by inertial-modes themselves ($' in eq. [I])- This so- 
called Cowling approximation is justified here. The po- 
tential perturbation 4>' is related to the density pertur- 
bation by the Poisson equation, 

V 2 $' = AxGp'. (70) 

Relative to the inertia term in equation {Tjl. the potential 
term is smaller by 



4nGp 



1 

X 2 ' 



(71) 



where we have used the results from equations @ and 
Q60(l . So the Cowling approximation is appropriate for 
high order modes (A > 1), which are indeed the modes 
we are concerned with (Paper II) 

We have assumed that turbulent viscosity does not 
modify mode structure significantly. This is equivalent 
of assuming that the turbulent dissipation rate 7 <C w. 
This assumption is confirmed by the numerical study to 
be presented in Paper II. 

For numerical tractability, we have ignored the rota- 
tional deformation of the planet, as well as the centrifugal 
force on the perturbed motion. These introduce ~ 10% 
error in the case of Jupiter, and much less for extra-solar 
planets. We do not expect the general structure and 
dispersion relation of inertial-modes to be significantly 
modified when these effects are taken into account. 

Our last simplifying assumption concerns the planet's 
atmosphere. This is by far our most uncertain as- 
sumption. In solving for the inertial-modes, we assume 
that the planet is fully convective (and neutrally buoy- 
ant). However, Jupiter- like planets have surface radia- 
tive zones with varying depths depending on their sur- 
face comp osition and ex t ernal irradiation. In the Jupiter 
model by Guillot ct al. (200J), the atmosphere is neu- 
trally stratified up to the photosphere (~ 1 bar in pres- 
sure), above which the temperature profile is largely 
isothermal with a scale height (~ depth) ~ 20 km ~ 
3 x 10 _4 R. Inertial-waves are reflected inward at a depth 
S ~ 1/A at most latitudes, and at a much shallower 
depth 6 ~ 1/A 2 near the singularity belt (JQJ. So 
modes with A < 1/3 x 10~ 4 ~ 3300 will have most of 
their WKB cavity inside the convective region and is only 
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mildly affected by the isothermal layer; while modes with 
A < 1/V3 x 10- 4 ~ 60 will have their upper-most WKB 
turning point below this isothermal layer and so will not 
be affected at all. In Paper II, we show that the modes 
relevant for tidal dissipation have A ~ 60. 

Two concerns arise following the above discussion. 

The referee (Dave Stevenson) called our attention to 
the problem of static stability in the Jovian atmosphere. 
This is currently uncertain due to the presence of wa- 
ter condensation. A 'wet adiabat' is less steep than a 
'dry adiabat' as water in an upward moving parcel con- 
denses and releases latent heat. So the actual tempera- 
ture profile will be super-adiabatic for some (very moist) 
parcels and sub-adiabatic for some (very dry) parcels. 
Such a 'conditional stability' makes it difficult to study 
the propagation of inertial-modes. 5 

But if the stably stratified atmosphere does extend suf- 
ficiently deep, and does harbor sufficiently strong buoy- 
ancy (Brunt- Vaisala N > u>), inertial-modes in the in- 
terior can couple to gravity-waves in the atmosphere to 
form hybrid modes. Since most of the mode inertia is in 
the interior, this will not affect much the frequencies of 
these modes; however, it may affect their surface struc- 
ture. Moreover, inertial gravity-waves in the atmosphere 
may propagate upward freely, break and dissipate at the 
low density environment in the upper atmosphere. This 
brings about an extra damping mechanism for the inte- 
rior inertial-modes. 

The second issue concerns extra-solar hot-jupiters. 
These planets are strongly irradiated by their host stars. 
Their surface isothermal layers may deepen to a depth 
of ~ 10 _2 i?. So even fairly low-order inertial-modes may 
be affected. 

For these two reasons, it is relevant to study behavior 
of inertial-modes in the presence of an isothermal layer. 
We plan to do so in the future. 

4.2. Special Case - R-modes 

R-modes have been considered as a potential candi- 
date for spinning down young neutron stars and for emit- 
ting d etectable gravitational waves (see, e.g. lOwen et alJ 
1998). So much attention have been paid to this special 
class of inertial-modes (they also appear in TablesQ&El). 
What is the relation betw een inertial-modes (also called 
'generalized R-modes' bv lLindblom fe Ipserlll99'9|) ') and 
R-modes? 

R-modes are purely toroidal, odd-parity, retrograde 
inertial-modes. They are, to the lowest order in 
Q, incompressible and mov e only on spherical shells 
(|Papaloizou fe Pringlel 11978(1 . Setting £ r = and V • 
£ = 0, we find that the displacement vector for R- 
modes can be expressed using a single stream function 
Q = Q(r, 9, <j>) where £ = V x (rQ), or 



Taking the curl of equation 1|11|) and retaining only terms 
to the lowest order in f£, we find 

™2 



sin 9 



Q, 



dQ 
' 89 



5 Data gathe red during the descent of the Galileo probe 
eiff et aT1ll99Sfl showed that, at the entry site, the temperature 
profile follows the dry adiabat from the photosphere (1 bar) down 
to 16 bars, with some evidence (temperature measurements and 
the presence of gravity waves) for a stable layer from 15 to 24 
bars. It is now known that the probe entered a particularly dry 
area, so this temperature profile may not be indicative of the whole 
atmosphere. 



1 8 



} 8Q 
89 



mq 



sm z 9 



Q = 0. (73) 



This is Legendre's equation with solution Q oc Ypj_ m (9, <p) 
and eigenvalue 11 = oj/2£1 = — m/£'(£' + 1) where £' > \m\ 
l|Papaloizou fe Pringlell9 784 . So the displacement vector 
£ = V x (rQ) oc r x VY^ iTn and is purely axial (toroidal) 
in nature. Obviously, both the angular dependence and 
the eigenfrequency of R-modes are independent of the 
equation of state (compare R-mode entries in Tablcs^fe 
EJ|. This is expected as R-modes are restricted to move 
only along spherical shells. 

According to this analysis, there should be infinite 
number of R-mode s for each m value . How ever, in ac- 
cordance with both lLindblom fe Ipserl |T999) and LF, we 
do not uncover any R- mode with £' > \m \. This is ex- 
plained by LF (see also lSchenk et al.ll2002l) . They real- 
ize that under the isentropic approximation, 6 poloidal- 
natured gravity-modes also have frequencies ~ 0(Q), 
besides from toroidal-natured R-modes. So these two 
branches of modes are allowed to mix and this elimi- 
nates the purely toroidal R-modes except for the lowest 
order one, £' = \m\. All other modes are a mixture of 
poloidal (made of terms depending on Yi^ m and VYp jm ) 
and toroidal (made of terms depending on xVY^' jm ) 
terms. These are our general inertial-modes. They can 
cause radial motion and their structure and their eigen- 
frequencies depend on the equation of state. 

Recall that we find eigenfunctions for uniform den- 
sity planets to be tp = ?/>i(xi)?/>2(£2) oc P™(xi) P™(x2)- 
When taking £ = \m\ + 1, we recover the above R-mode 
solution. 

Since R-modes are insensitive to the equation 
of state, they can exist even in radiative stars 
ijPapaloizou fc Pringlel 119781) . In contrast, general 
inertial-modes exist only in neutrally buoyant medium. 

5. SUMMARY 

In this work, we have studied inertial-modes with the 
purpose of unraveling the role they may play in the tidal 
dissipation process of Jupiter. 

With the hel p of the ellips oidal coordinates (0:1,352) 
first adopted bv iBrvanl l)1889|) . we have shown that the 
partial differential equation governing inertial-modes in 
a sphere of neutrally buoyant fluid can be separated into 
two ordinary differential equations when the density is 
uniform or when the density has a power-law dependence. 
The latter case is a novel result as far as we know. For 
more general density scalings, we show that we can ob- 
tain an approximate solution to the inertial-modes that 
is accurate to the second order in wave- vector. This im- 
portant result underlies much of our analytical study. 

The dispersion relation /1 = uj/2fl w sin(ri27r/2(rii + 
n 2 )) (eq. |3£] ) rather generally describes how the fre- 
quency of an inertial-mode depends on its structure. The 
quantum numbers (711,712) are respectively the number 
of nodes in the x\ and X2 coordinates, and the third 
quantum number m is the conventional azimuthal num- 
ber. In our notation, positive m denotes prograde modes, 

6 Isentropic refers to the fact that both the background model 
and its adiabatic perturbation satisfy the same equation of state. 
For instance, adiabatic perturbations in neutrally buoyant models. 
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and negative m retrograde modes. So frequencies of 
inertial-modes depend on the direction of wave propaga- 
tion, with modes propagating close to the rotation axis 
having higher frequencies. This dispersion relation also 
indicates that inertial-modes are dense - for any given 
frequency, one can always find a combination of n\ & ni 
that approaches it sufficiently closely. 

We find that inertial-modes naturally cause small 
(but non-zero) Eulerian density perturbations, of order 
wTdyn/A (where Td yn is the dynamical time-scale) smaller 
than those by p-modes of comparable displacement am- 
plitude (eq. jSSl)- Their motion is nearly anelastic. This 
implies that inertial-modes are only weakly coupled to 
the tidal potential. It also implies that inertial-modes 
are not dissipated primarily through heat diffusion, but 
rather through viscosity. 

In its propagating region, an inertial-mode has a wave- 
vector k (measured in x%,X2 coordinates) that is nearly 
constant and is insensitive to the local scale height and 
density distribution (eq. |46j ) . In this region, the ampli- 
tude envelope of an inertial-mode rises with lowering of 
the density as An inertial-mode encounters its up- 

per turning point when the density scale-height becomes 
small and comparable to the wave-vector. The depth of 
this turning point depends on the latitude: when away 



from the spherical angle 6 = cos^ 1 ±/x, the turning point 
occurs at a depth ~ R/X; while near this angle, reflec- 
tion occurs at a much shallower depth ~ R/X 2 . Here, 
A ~ 2(ni+n2)- We call the special surface region around 
| cosf?| = /j, the "singularity belt". An inertial-mode has 
the highest amplitude as well as the sharpest spatial gra- 
dient inside this belt. This region is associated with the 
strongest turbulent dissipation. 

Among the many simplifying assumptions we have 
adopted, the most uncertain one concerns the static sta- 
bility in the planet atmosphere ( iJ4.1[l . We plan to study 
its effects on inertial-mode structure and dissipation in a 
future work. 
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APPENDIX 
ELLIPSOIDAL COORDINATES 

iBrvanl <)1889fl introduced a new set of coordinates under which the left-hand-side of equation (fT4")l becomes separable. 
For any given density profile, the right-hand-side of the same equation is generally inseparable - however, for the two 
special cases of uniform and power-law density profiles ( H2.2I and !2.3fl . the right-hand-side becomes separable and we 
can obtain (semi)-analytical solutions for the inertial-modes. This new set of coordinates allow us to obtain good 
approximate solutions for more general density profiles ( H2.4|l . 

We scale all length in the problem by the radius of the planet so that the Cartesian coordinates, x, y and z, fall 
within the range of (—1, +1). The new ellipsoidal coordinates (x\, X2, </>) are related to the Cartesian coordinates as 




XlX 2 



with xi S [/x, 1], X2 € 
are given by 



-/x, /j], and (p is the usual azimuthal angle with 

(l-xl)(l-xl) 



(Al) 

£ [0,27r]. The cylindrical and spherical radii 



y = 



(i-M 2 ) 



Wu 



TABLE 1 

INERTIAL-MODE EIGENFREQUENCIES a IN A UNIFORM DENSITY 
SPHERE. 



{(. - \m\) parity 15 |m| = 1 |m| = 2 |m| = 3 |m| = 4 



l c 


o 


1 nnnn 


-U.UUU 1 


-u.ouuu 




2 


e 


-1.5099 


-1.2319 


-1.0532 


-0.9279 




e 


0.1766 


0.2319 


0.2532 


0.2613 


3 


o 


-1.7080 


-1.4964 


-1.3402 


-1.2203 




o 


-0.6120 


-0.4669 


-0.3779 


-0.3175 




o 


0.8200 


0.7633 


0.7181 


0.6807 


4 


e 


-1.8060 


-1.6434 


-1.5119 


-1.4044 




e 


-1.0456 


-0.8842 


-0.7735 


-0.6920 




e 


0.0682 


0.1018 


0.1204 


0.1312 




e 


1.1834 


1.0926 


1.0222 


0.9652 


5 


o 


-1.8617 


-1.7340 


-1.6236 


-1.5290 




o 


-1.3061 


-1.1530 


-1.0401 


-0.9525 




o 


-0.4404 


-0.3595 


-0.3040 


-0.2635 




o 


0.5373 


0.5100 


0.4869 


0.4669 




o 


1.4042 


1.3080 


1.2309 


1.1670 



a Here, we present cigenfrequencies as 2/im/|m| = 
SIGN[m] cj/f2, where u> is the eigenfrequency in the rotating 
frame, and Q the rotational frequency. So positive values de- 
note prograde modes, while negative ones retrograde (opposite 
to that in LF). The eigenvalues t and m appear in the eigen- 
function as ip = P™ (x\) P™ (xz) ■ For each (£, |m|) pair, there 
are (£— \m\) distinct eigenfrequcncies. These results agree with 
those obtained by LF. 

b This denotes the parity with respect to the equatorial plane: 
e for even and o for odd. 

c This row shows the eigenfrequcncies of pure r-modes. They 
are odd-parity, retrograde modes satisfying ui = 2Q/(\m\ + 1) 
to the first order in Q. 



TABLE 2 

ElGENFREQUENCIES a OF \m\ = 2 INER.TIAL-MODES FOR VARIOUS POWER-LAW DENSITY 

PROFILES 



(<-M) 


parity 


13 = 0.0 


(3 = 0.1 


(3 = 0.5 


(3 = 1.0 


(3 = 2.0 


p = fcp 2 (LF) b 


i c 


o 


-0.6667 


-0.6667 


-0.6667 


-0.6667 


-0.6667 


-0.6667 


2 


e 


-1.2319 


-1.2133 


-1.1607 


-1.1224 


-0.8628 


-1.1000 




e 


0.2319 


0.2673 


0.3830 


0.4860 


0.6159 


0.5566 


3 


o 


-1.4964 


-1.4789 


-1.4256 


-1.3822 


-1.3317 


-1.3578 




o 


-0.4669 


-0.4732 


-0.4924 


-0.5082 


-0.5270 


-0.5173 




o 


0.7633 


0.7942 


0.8919 


0.9761 


1.0798 


1.0259 


4 


e 


-1.6434 


-1.6287 


-1.5820 


-1.5415 


-1.4909 


-1.5196 




e 


-0.8842 


-0.8811 


-0.8726 


-0.8671 


-0.8628 


-0.8629 




e 


0.1018 


0.1188 


0.1780 


0.2364 


0.3199 


0.2753 




e 


1.0926 


1.1137 


1.1814 


1.2408 


1.3150 


1.2729 


5 


o 


-1.7340 


-1.7217 


-1.6820 


-1.6460 


-1.5988 


-1.6272 




o 


-1.1530 


-1.1471 


-1.1286 


-1.1133 


-1.0954 


-1.1044 




o 


-0.3595 


-0.3669 


-0.3904 


-0.4108 


-0.4367 


-0.4217 




o 


0.5100 


0.5313 


0.6020 


0.6673 


0.7554 


0.7039 




o 


1.3080 


1.3226 


1.3700 


1.4122 


1.4654 


1.4339 



a Eigenfrcqucncics arc again presented as 2/im/|ra| = SIGN[ra] with positive 

values denoting prograde modes, and negative ones retrograde modes. Here we define 
I = 2(m + U2) + \fn\ — 5, where <5 = for even-parity and 1 for odd-parity, and ni, 
ri2 are the number of nodes in x\ and X2 ranges, respectively. When (3 = 0, this t 
matches that appearing in the Legendre polynomial (-PI™). The number of eigenmodes 
is conserved when f3 varies, with a close one-to-one correspondence between modes in 
different density profiles. 

b For comparison, we list results calculated by LF (their Table 5) using a series expan- 
sion method for a polytrope model p <x p . Not surprisingly, their eigenfrequcncies fall 
somewhere in between our results for (3 = 1 and (3 = 2 models. 

c This row shows pure r-modes. Both the frequency and the eigenfunction of these 
modes do not depend on the equation of state, to the lowest order in Q. 
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Fig. 1. — The ellipsoidal coordinates for p = 0.75. Left panel: equidistant curves of constant x\ (solid curves) and constant X2 (dashed 
curves) in a meridional plane. They become dotted curves outside the surface of the planet, which corresponds to either xi = p or | rca ] = fjb 
or both (at the special latitude | cos0| = p). The equidistant curves are most closely packed near the surface at this special latitude. Over 
the rest of the sphere, the x\ axis runs similarly as the cylindrical radius, while the X2 axis runs largely parallel to the rotation axis. Right 
panel: equidistant curves of constant spherical radius in the x\ — X2 plane. 




Fig. 2. — Left-hand plot shows the logarithm of the density (logp) as a function of the logarithm of the depth (log(l — r)) for a /3 = 1 
power-law model (solid curve). The density scale is arbitrary. The dotted line is the density profile for a n = 1 polytrope (p <x p 2 ). The 
two profiles behave similarly except near the center. Right-hand panel: eigenfunctions an d ^2(22) as a function of x± or X2 for an 

even-parity, retrograde m = — 2, p = 0.4881 mode in the f3 = 1 model (the solid curve). The global eigenfunction ip(xi,X2) = ipi ip2- The 
dashed line demarcates the x\ and X2 boundary at p = 0.4881. There are 5 and 3 nodes within x\ and X2 ranges, respectively, with a 
corresponding i = 2(ni + 122) + \m\ = 18 (see Table l2l . As Xi approaches p, the solid curve exhibits a rise in amplitude that is described 
by a 1/^/p WKB-envelope. For comparison, we draw (in dotted curve) the eigenfunction of a counterpart mode in the /3 = (uniform 
density) model. 



2 1 2 1 2 

: x + y + z 



(1 - M >2 ■ 

The Jacobian that relates the new coordinate system (xi,X2,4 > ) to the Cartesian (x,y,z) is 



J = det 

so the volume element dx dy dz = Jdx\ dxi 



dx/dxi dx/dx2 dx/dcj) 
dy/dxi dy/dx2 dy/d<p 
dz/dxi dz/dx2 dz/dcf) 
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Fig. 3. — Mode energy as a function of ri\ for m = 
are normalized with ip2{x2 = 0) = 1 and satisfy n\ 



—2 even-parity inertial-modes in a (3 = 1 model (solid curve). Here, all inertial-modes 
712 (so p ~ 0.70). The dotted curve is a fitting power-law E <x n| ,6s . 




Fig. 4. — A meridional look at the inertial-mode presented in Fig. |3 This mode is symmetric with respect to the equator, retrograde 
(m = —2), with ni = 5, 712 = 3. The planet model has a density profile /3 = 1. The left panel shows the Eulerian density perturbation 
(p'), using both gray-scale and equidistant contours. Lighter regions (and dashed contours) stand for p' < 0, while darker regions (solid 
contours) represent p' > 0. The dotted curves indicate the (21,2:2) coordinates. Counting the number of nodes along each coordinate, one 
recovers n\ and 712. The right panel shows the fluid velocity (v r & vg components only) in the rotating frame as ar rows , with the size of 
the arrows proportional to ^/v. Notice that v r vanishes at the surface, as is required by the boundary condition (eq. 1271 '). Both the mode 
amplitude and the wave-vector remain relatively constant over much of the planet, but rise sharply toward the surface. This rise is most 
striking near the special angle | cos£>| = p. = 0.4881, marked here by straight lines. 



Fig. ^depicts the equi-distance curves of (x±, X2) in a meridional plane, as well as the equi-distance curves of radius r 
in the X\ —x-i plane. In 3-D, surfaces of constant x\ appear as co-axial cylinders around the z-axis (or prolate ellipsoids), 
while surfaces of constant Xi resemble bandannas symmetrical with respect to the equator (or oblate ellipsoids). On 
the spherical surface, either x\ or 1 0:2 1 (or both) equals /j, - we call the region where x± « \x2[ ~ [i the 'singularity 
belt', a region of special significance. Moreover, x\ = 1 at the rotation axis and X2 = at the equator. 

Partial differentiation with respect to vj, z and r can be expressed in the new coordinates as, 
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Fig. 5. — Surface density perturbation (p') for the same mode (m = — 2,ni = 5,rt2 = 3) as in Fig. [I] The short vertical stick marks the 
rotation axis. In both hemispheres, there is a latitudinal belt within which density perturbation attains the largest amplitude and varies 
with the steepest gradient. This lies around | cos0| = /i and we call it the "singularity belt". The surface pattern of the mode also discloses 
the quantum numbers: from the pole to the equator, there are n\ nodal patches above the singularity belt, and n,2 nodes below it. 




Fig. 6. — As Fig. |1] except for a m = — 2 R-mode. This mode has fi = l/(\m\ + 1) = 0.3333 and in our notation, n\ = ri2 = and 
£ = \m\ + 1 = 3. There is no internal nodal point and the radial velocity is zero everywhere. This mode exhibits odd symmetry toward the 
equator. 
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Fig. 7. — Frequencies of even-parity, \m\ = 2 inertial-modes (shown here as SIGN/i) as a function of wave-number A ~ 2(n\ + n,2). Modes 
above the dotted line are prograde modes, and retrograde if below. These frequencies are obtained using a uniform density model but the 
overall features persist for mo re g eneral models. Note the asymmetry between the frequencies of prograde and retrograde modes — this 
follows from equations 1361 & 1371 . 



MAKING (3 MODELS 

We construct a hydrostatic model for a power-law density profile of the form p = (1 — r 2 ) 13 (eq. [HH])- Any value of 
P is allowed, in contrast to the polytrope case where the polytrope index has to be smaller than 5. 

For our numerical integration, we adopt lnp as the independent variable, and hip and M(r) as the variables. Here 
p is pressure and M(r) is the mass within radius r. The normal hydrostatic equations apply with the boundary 
conditions that at the center M(r) = and at the surface, 



P 



surface 



(l-r) 



/3+1 



(Bl) 



where surface gravity gsurfacc = GM(r = l)/r(r = l) 2 . This is obtained using dp/dr = —pg. We can then solve the 
boundary value problem to obtain the interior structure of such a model. Here, r runs from to 1 and the density 
scale is arbitrary. 



